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I.  Introduction 


A  weakness  of  many  operational  numerical  weather  prediction  (NWP)  mod¬ 
els  has  been  the  forecast  accuracy  of  the  distribution  and  intensity  of  pre¬ 
cipitation,  and  the  accompanying  forecast  fields  of  atmospheric  moisture. 
Most  models  use  the  equation  for  the  conservation  of  mass  of  water  vapor 
represented  by  some  form  of  absolute  humidity,  along  with  source  and  sink 
parameterizations  of  evaporation-condensation  processes  such  as  large-scale 
and  convective  scale  precipitation.  Certainly,  the  parameterizations  may 
inadequately  simulate  the  corresponding  atmospheric  phenomena,  hut  improper 
specification  of  the  initial  distribution  of  atmospheric  moisture  could  also 
contribute  significantly  to  forecast  errors.  Many  presently  used  moisture 
analysis  methods  treat  moisture  in  much  the  same  manner  as  the  other  mass 
fields,  without  any  special  attention  to  any  similarities  or  differences 
that  may  exist  between  moisture  and  other  mass  field  distributions.  This 
report  describes  a  study  which  attempts  to  determine  the  general  character¬ 
istics  of  the  global  moisture  distribution  and  to  compare  the  features  of 
their  spectral  composition  to  those  of  the  mass  and  motion  field  parameters. 
This  comparison  will  allow  for  the  more  realistic  design  of  a  moisture  anal¬ 
ysis  method  in  which  the  similarities  and  differences  of  moisture  with  mass 
and  motion  fields  can  be  accounted  for.  The  result  should  be  a  more  accu¬ 
rate  representation  of  the  initial  moisture  field. 

Controversy  has  arisen  over  the  importance  of  the  initial  specifica¬ 
tion  of  moisture  in  precipitation  and  moisture  distribution  forecast  accura¬ 
cy.  In  a  study  by  Phillips  (1978)  \  two  fine  mesh  limited  area  models  were 
run  with  two  different  moisture  analyses,  one  of  which  fit  the  radiosonde 
values  of  moisture  much  better  in  an  area  adjacent  to  the  Gulf  of  Mexico 
than  the  other.  In  both  models  much  more  realistic  forecasts  of  vorticity 
and  precipitation  resulted  in  the  Gulf  area  from  using  the  analysis  which 
was  more  representative  of  the  observations.  Apparently,  the  specification 
of  too  large  values  for  humidity  at  mid-tropospheric  levels  resulted  in  the 
overdevelopment  of  cloudiness  which  in  conjunction  with  the  associated  larger 
release  of  latent  heat  causad  an  overdevelopment  of  vorticity  and  excess 

1.  Phillips,  R,  A.,  1978«  A  Test  of  Finer  Resolution.  Office  Note  171, 

National  Meteorological  Center,  National  Weather  Service,  NOAA. 
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precipitation,  van  Maanen  (1979)  describes  an  experiment  conducted  at  the 
Europeon  Centre  for  Medium  Range  Forecasts  (ECMWF)  in  which  three  different 
sets  of  initial  conditions  are  used  in  a  five-day  forecast  using  the  ECMWF 
spectral  model.  The  three  sets  of  conditions  were  (a)  mass,  wind  and  humid¬ 
ity  fields  from  the  National  Meteorological  Center  (NMC)  analysis,  (b)  the 
ECMWF  analysis  for  the  same  fields,  and  (c)  the  ECMWF  in  which  the  relative 
humidity  was  set  at  60  percent  for  all  model  levels.  Initial  state  (c)  re¬ 
sulted  in  differences  from  (b)  in  five-day  forecast  fields  of  pressure  and 
humidity  that  were  of  the  same  size  as  differences  between  forecasts  using 
(a)  and  (b) .  The  experiments  were  conducted  for  only  one  date,  but  assuming 
that  the  results  hold  for  other  dates,  it  was  concluded  that  the  specifica¬ 
tion  of  initial  moisture  is  important  out  to  at  least  five  days  of  forecast 
time.  By  changing  from  an  omni-directional  successive  correction  humidity 
analysis  to  one  which  took  into  account  the  direction  of  humidity  gradients, 
Atkins  (1974) 3  demonstrated  significant  improvement  in  rainfall  forecasts 
from  a  10-level  model  for  a  specific  weather  event  over  western  Europe.  Her 
conclusion  was  that  the  quality  of  numerical  short  range  forecasts  of  precip¬ 
itation  depends  highly  on  the  correctness  of  the  initial  specification  of  the 
humidity  fields. 

Somewhat  different  conclusions  have  been  reached  as  a  result  of  other 

4 

studies.  Smagorinsky  et  al.  (1970)  attempted  to  determine  which  atmospher¬ 
ic  variables  had  to  be  specified  initially  and  which  could  be  derived  either 
diagnostically  or  were  secondary  in  importance  because  of  dynamic  coupling 
in  forecast  models.  Using  the  Geophysical  Fluid  Dynamics  Laboratory  (GFDL) 
nine-level  hemispheric  model  for  a  particular  set  of  initial  conditions,  a 
two-week  integration  was  carried  out.  Then  to  test  the  importance  of  speci¬ 
fication  of  initial  humidity,  another  two-week  integration  was  performed  in 
which  the  only  change  in  the  initial  conditions  was  to  assign  the  relative 

2.  van  Maanen,  J.,  1979s  Description  of  Some  Properties  of  the  Humidity 

analysis  Scheme  in  Use  at  the  European  Centre  for  Medium  Range  Weather 
Forecasts.  Verslagen,  V-334,  Koninkluk  Nederlands  Meteorologisch  In- 
stituut.  (Available  from  NTIS) 

3.  Atkins,  M.  J.,  1974s  The  objective  analysis  of  relative  humidity. 

Tellus,  26,  663-671. 

4.  Smagorinsky,  J.,  K.  Miyakoda,  and  R.  F.  Strickler,  1970:  The  relative 

importance  of  variables  in  initial  conditions  for  dynamical  weather 
prediction.  Tellus,  22,  141-157. 
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humidity  a  value  of  70  percent  everywhere.  Their  results  indicated  that  in 
extra-tropical  latitudes  there  were  no  significant  differences  in  the  large 
scale  variations  of  both  relative  humidity  and  accumulatad  precipitation 
after  about  12  hours  of  forecast  time.  Their  conclusion  wa6  that  a  proper 
specification  of  the  horizontal  wind  field  in  the  extra- tropical  latitudes 
will  set  up  those  vertical  velocity  patterns  that  play  the  major  role  in  the 
spatial  distribution  of  moisture  and  precipitation  in  forecasts  exceeding 
about  12  hours.  Lejenas  (1979) 5,  in  experiments  with  a  primitive  equation 
model  covering  the  Northern  Hemisphere  down  to  30°N,  concluded  that  moisture 
initialization  was  only  really  important  in  assuring  precipitation  forecast 
quality  in  the  first  six  hours  of  the  forecast.  His  results  suggested  that 
vertical  velocities  implicit  within  the  initial  horizontal  wind  field  are  at 
least  as  important  as  fine  detail  in  the  moisture  fields,  especially  for 
short  term  (<24  hour)  forecasts,  and  that  moisture  specification  is  of  little 
importance  in  longer  term  (>48  hour)  forecasts  for  extra 'tropical  regions. 

No  attempt  is  made  in  this  study  to  resolve  the  question  of  the  impor¬ 
tance  of  initial  specification  of  moisture  in  numerical  forecast  models. 
Perhaps  part  of  the  apparent  discrepancies  in  the  studies  mentioned  above 
could  be  attributed  to  the  individual  cases  studied,  or  more  likely,  to  the 
differences  in  the  forecast  models  used.  An  even  more  important  problem, 
however,  is  the  consistent  deficiencies  seen  in  precipitation  forecasts  from 
operational  models  which  perform  well  in  forecasting  mass  and  motion  fields 
over  several  days.  Why  dynamic  coupling  does  not  result  in  a  consistently 
accurate  forecast  of  precipitation  in  these  models  is  a  question  that  has 
not  been  satisfactorily  answered.  From  the  diversity  of  results  from  pre¬ 
vious  studies  mentioned  above,  one  may  conclude  that  initial  moisture  spec¬ 
ification  is  more  important  in  seme  models  than  it  is  in  others.  Only  tests 
performed  in  conjunction  with  a  particular  model  over  many  cases  can  lead  to 
concrete  conclusions  about  the  importance  of  initial  moisture  specification 
in  that  model. 

If  the  initial  specification  of  the  moisture  field  is  important  for 
quality  humidity  and  precipitation  forecasts  from  at  least  some  models,  it 
seems  that  the  most  generally  desirable  objective  would  be  to  fashion  the 
initial  moisture  field  to  most  closely  resemble  the  highest  quality  moisture 
observations.  The  mark  of  a  good  analysis  scheme  for  a  given  model  is  the 

5.  Lejenas,  H.,  1979:  Initialization  of  moisture  in  primitive  equation 

models.  Mon.  Wea.  Rev.,  107,  1299-1305. 


consistent  production  of  the  best  forecast  that  can  be  achieved  by  that  mod¬ 
el.  A  good  representation  of  the  observations  is  an  important  requirement 
that  an  analysis  scheme  should  meet  in  preparing  data  for  the  models  that 
best  simulate  atmospheric  processes.  Model  specific  modifications  may  be 
necessary  after  this  point,  but  the  necessity  of  such  modifications  would  be 
more  a  reflection  on  the  reliability  of  the  model  than  on  that  of  the  analy¬ 
sis.  In  fact,  it  may  be  more  advantageous  to  modify  or  redesign  a  moisture 
analysis  scheme  once  the  properties  of  the  spatial  distribution  of  moisture 
are  known,  so  that  the  scheme  could  consistently  produce  an  analysis  which 
reflects  the  unique  variations  within  the  moisture  field  that  may  prove  im¬ 
portant  in  preserving  the  accuracy  of  the  predicted  moisture  and  precipita¬ 
tion  field. 

For  a  global  spectral  forecast  model,  the  proper  specification  of  the 
wave  structure  of  the  moisture  field  is  likely  to  be  of  greater  importance 
than  point-wise  agreement  with  the  observations.  At  the  Air  Force  Geophysics 
Laboratory,  a  global  spectral  model  of  the  moist  atmosphere  is  being  devel¬ 
oped.  In  this  model,  particular  importance  will  be  placed  on  moisture  phys¬ 
ics  and  initialization,  in  an  attempt  to  produce  more  skillful  cloud  fore¬ 
casts  in  support  of  Ai*  Force  cloud-free  line  of  sight  requirements.  There¬ 
fore,  proper  specification  of  initial  moisture  distribution  will  be  an  im¬ 
portant  prerequisite  for  spectral  model  integrations.  To  know  how  an  ob¬ 
jective  analysis  scheme  must  be  designed  to  best  represent  the  unique  spec¬ 
tral  characteristics  of  the  spatial  moisture  variation,  the  spectral  charac¬ 
teristics  of  the  global  moisture  distribution  must  be  known.  This  study  ex¬ 
amines  the  zonal  and  spherical  harmonics  of  mass,  motion,  and  moisture  fields 
over  periods  of  one  week  of  summer  data  and  one  week  of  winter  data.  The  re¬ 
sults  obtained  from  the  study  may  then  be  applied  in  designing  an  analysis 
procedure  that  takes  into  account  both  the  conmon  and  unique  properties  of 
global  moisture  distributions. 

II.  Description  of  Study 

In  a  previous  study  (Gerlach,  19816) ,  two  synoptic  times  were  chosen 
for  spectral  analysis.  The  present  study  expands  the  data  sample  to  two  one- 
week  periods*  January  15  (OOGKT)  through  January  21  (12GMT) ,  1978  (referred 

6.  Gerlach,  A.  M.,  Ed.,  1981s  Technique  Development  for  Weather  Forecasting. 

AFGL-TR-82-0020,  Final  Report,  Contract  F19628-81-C-0039,  Systems  and 

Applied  Sciences  Corporation  (SASC),  pp.  8-92. 
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to  as  JAN),  and  July  16  (OOGMT)  through  July  22  (12GMT),  1978  (JUL) .  Global 
observations  and  analyses  for  OOGMT  and  12GMT  for  each  day  in  both  periods 
were  extracted  from  data  tapes  from  the  First  GARP  Global  Experiment  (FGGE) 
ot  the  Global  Atmospheric  Research  Program  (GARP) .  Level  II-A  and  Level 
1II-A  data  sets  were  used  for  the  observations  and  the  2.5°  latitude- longi¬ 
tude  gridded  analyses,  respectively.  These  "A"  sets  contain  only  those  ob¬ 
servations  that  were  received  within  a  preset  cutoff  time  and  used  for  real¬ 
time  FGGE  purposes.  In  the  present  study,  the  Level  II  observations  were 
used  to  "correct"  the  Level  III  analysed  fields  in  order  to  produce  a  grid 
point  representation  of  the  observations  at  each  synoptic  time.  The  method 
and  motivation  for  the  corrections  will  be  described  subsequently. 

For  each  synoptic  time,  only  radiosonde  and  pilot  balloon  observations 
from  the  Level  II-A  tapes  were  used.  In  addition,  only  those  observations 
were  used  that  satisfied  the  quality  control  checks  performed  on  the  FGGE 
data  before  they  were  put  on  tape.  The  resulting  FGGE  Level  II-A  observa¬ 
tion  values  were  plotted  on  mercator  projection  maps  for  1/15/78  and  7/20/78, 
both  at  OOGMT.  Examples  of  these  maps  including  contours  from  hand-analysis 
are  shown  in  Figs,  la-lf.  The  maps  are  presented  only  to  show  the  key  synop¬ 
tic  features  in  the  Northern  Hemisphere  for  those  times.  The  contours  drawn 
over  oceans  and  in  the  Southern  Hemisphere  are  oniy  estimates  due  to  the 
sparsity  of  observations  in  these  regions.  The  synoptic  cases  shown  are 
typical  of  circulations  observed  in  the  two  seasons,  and  as  such  are  repre¬ 
sentative  of  normal  seasonal  conditions... 

Before  considering  specific  characteristics  of  these  data  sets,  we 
first  review  several  well  known  properties  of  global  moiBture  fields.  The 
zonally  averaged  structure  of  the  specific  humidity  (q) ,  for  example, 
always  reflects  that  of  climatology:  that  is,  zonally  averaged  q  (q) 
is  a  maximum  at  the  equator  and  decreases  monotonically  toward  both  poles. 
Departures  from  the  zonal  mean  for  any  given  observation  result  primarily 
from  either  the  longitudinal  land  mass  distribution  or  disturbances  in  the 
zonal  mean  flow.  The  magnitude  of  the  departures  due  to  ocean-land  distri¬ 
bution  can  be  obtained  from  climatological  studies. 

g 

Peixoto  (1970)  conducted  such  a  climatological  study  in  which  he  plot- 

7.  Obtained  from  Department  of  the  Air  Force,  OL-A,  USAF  Environmental  Tech¬ 

nical  Applications  Center  (MAC),  Federal  Building,  Asheville,  NC  28801. 

8.  Peixoto,  J.  P.,  1970:  Hater  vapor  balance  of  the  a tgp sphere  from  five 

years  of  hemispheric  data.  Hordic  Hydrol.,  2j  120-136. 
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Fig.  la.  Subjective  analysis  of  mean  sea  level  pressure  (mb)  for  1/15/78  at 
OOGMT.  Only  rawinsonde  and  surface  land/marine  observations  were  used  and 
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ted  the  time-averaged  vertically  integrated  precipitable  water  defined  in 
terms  of  the  time-averaged  specific  humidity  §  as 

300 

ft  =>  g  f  q  dp  (1) 

*sfc 

where  g  is  the  acceleration  due  to  gravity  and  p  is  pressure.  His  plot, 

9 

taken  from  Rosen  et  al.  (1979)  ,  is  reproduced  in  Fig.  2.  By  removing  the 
zonal  mean  at  each  latitude  from  both  the  time-averaged  field  and  the  synop¬ 
tic  field,  the  departure  due  to  orientation  of  ocean-land  masses  can  be  com¬ 
pared  with  that  due  to  eddies.  When  this  was  done  for  the  1/15/78  OOGMT 
case,  calculations  using  q  data  for  a  station  at  43.72°N  and  65.25°W  (near 
the  moisture  maximum  associated  with  the  northwest  Atlantic  storm)  indicated 

a  departure  of  precipitable  water  from  the  zonal  mean  for  45°N  (computed 

-2 

from  corrected  q  fields  to  be  described  later)  of  2.3  x  10  kg  m  .  The 
largest  departure  from  the  zonal  mean  of  the  time-averaged  values  on  Fig.  2 
occurs  at  about  30°N  (over  the  Sahara)  and  was  calculated  at  -1.8  x  10  kg  m  2 
(using  a  time-averaged  zonal  mean  computed  from  Oort  and  Rasmusson  (1971) 10 
for  30°N,  yearly  average) .  From  this  comparison  it  is  clear  that  eddy  de¬ 
partures  can  be  as  large  as  the  largest  time-averaged  departures  resulting 
from  ocean-land  distribution. 

The  moisture  variables  examined  in  this  study  were  specific  humidity 
and  layer  precipitable  water.  Relative  humidity  was  not  used  because  of  its 
much  higher  degree  of  spatial  irregularity.  Initially,  maps  such  as  those 
in  Fig.  1  were  prepared  for  relative  humidity,  and  it  was  obvious  that  the 
observations  were  characterized  by  many  pockets  of  high  humidity  located  ad¬ 
jacent  to  areas  of  low  humidity.  Such  a  discontinuous  field  is  very  diffi¬ 
cult  to  analyze  subjectively,  and  the  ability  of  any  objective  analysis 
scheme  to  replicate  such  a  highly  irregular  field  is  questionable. 

Radiosonde  observations  of  both  temperature  T  and  dewpoint  depression 
0^  at  mandatory  and  significant  pressure  levels  were  used  in  this  study  to 
calculate  Level  II  values  of  q.  To  obtain  corresponding  values  of  W,  the 
calculated  values  of  q  up  to  250  mb  were  used  in  a  numerical  integration  ap¬ 
proximating  the  vertical  integrals 

9.  Rosen,  R.  D.,  D.  A.  Salstein,  and  J.  P.  Peixoto,  1979:  Variability  in 

annual  fields  of  large-scale  atmospheric  water  vapor  transport.  Hon. 
Wea.  Rev.,  107,  26-37. 

10.  Oort,  A.  H.,  and  E.  M.  Rasmusson,  1971:  Atmospheric  Circulation  Sta¬ 

tistics.  NOAA  Prof.  Pap.  No.  5  (NTIS  Com-72-50295) ,  323  pp. 
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Fig.  2.  Polar  stereographic  map  of  vertically  integrated  time  mean 
precipitate  water  W  for  the  period  Hay  1,  1958  -  April  30,  1963. 
Units  for  w  are  x  10  kg  m”  (from  Rosen  et  al.,  1979  »  taken  origin¬ 
ally  from  Peixoto,  1970  ; . 


“sso  -  5^  /  q  dp 

sfc  -  300 

for  each  sounding  (p.  is  the  density  of  water,  and  W_Crt  and  are  computed 

U  850  500 

in  cm) .  The  subscripts  850  and  500  are  used  only  to  indicate  "lower"  and 
"upper"  layers;  in  most  soundings  they  represent  approximately  the  mid-levels 
of  the  two  layers.  The  layers  were  chosen  so  that  the  mid-levels  would  coin¬ 
cide  with  the  pressure  levels  at  which  the  harmonic  analysis  was  performed 
on  the  other  quantities  considered  in  this  study. 

To  insure  a  sufficiently  complete  set  of  q  values  for  the  above  integra¬ 
tion  at  each  sounding  location,  the  following  procedure  was  used.  First,  a 
check  was  made  to  determine  if  a  q  value  existed  at  a  pressure  less  than  or 
equal  to  300  mb.  If  not,  a  check  was  made  to  see  if  the  highest  level  with 
a  q  value  in  that  sounding  was  between  300  and  400  mb.  If  so,  a  value  was 
extrapolated  to  300  mb  according  to  the  formula 

q300  =  q0  Cl  "  *695  <p0  "  300>/100:j  O) 

where  qQ  is  the  specific  humidity  value  at  the  highest  pressure  level  pQ. 

Using  specific  humidity  values  at  300  and  400  mb  from  Oort  and  Rasmusson 
(1971) representing  zonal,  yearly  averages,  the  quotients  q3Q0  /  q4Q0  were 
averaged  over  5°  latitude  intervals  from  0°  to  75°N.  An  average  value  of 
.305  was  obtained,  from  which  the  complement  value  of  .695  was  obtained.  A 
linear  variation  was  assumed  between  the  two  levels.  At  this  point  a  test 
was  performed  to  check  if  (a)  the  lowest  pressure  was  sPgfc  ~  300  mb,  or  (b) 
at  least  one  value  had  a  pressure  £  300  mb  (including  the  extrapolated  value 
if  one  had  been  added) ,  or  (c)  both.  The  sounding  was  discarded  if  none  of 
the  above  was  true.  For  type  (a),  there  had  to  be  at  least  three  values 

>pgfc~  300  mb,  and  the  lowest  p  value  had  to  be  £  700  mb  in  order  for  the 

sounding  to  be  used.  Case  (b)  required  at  least  three  values  with  pressures 
such  that  300  £p£  Pafc~  300  for  the  sounding  to  be  used.  A  minimum  of 
three  levels  was  also  required  for  case  (c) .  A  cubic  spline  was  then  used 

to  interpolate  values  of  q  at  even  multiples  of  50  mb  between  the  highest  and 

lowest  pressure  value  in  the  sounding.  Using  the  given  and  interpolated  val¬ 
ues  (and  extrapolated  300  mb  values  where  they  exist),  the  integrals  (Eq.  (2)) 
were  evaluated  using  the  trapezoidal  rule.  When  negative  values  of  q  were  gen¬ 
erated  from  the  cubic  spline  (which  generally  occurred  in  less  than  10  percent 


W 


500  gp 


sfc 


300 


300 
q  dp 


(2) 
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of  the:  accepted  soundings  for  an  observation  period  in  this  study),  the  nega¬ 
tive  values  were  replaced  by  assuming  a  linear  variation  of  q  with  pressure 
between  the  non-negative  values  immediately  above  and  below  the  negative 
value  or  values. 

The  same  basic  scheme  was  used  to  calculate  W„, „  and  for  Level  III 

850  500 

values,  except  in  this  case  a)  analyzed  values  of  T  and  RH  were  available 
to  calculate  q  at  mandatory  pressure  levels  only,  and  b)  the  surface  pres¬ 
sure  values  needed  in  Eq.  (2)  had  to  be  derived.  To  derive  surface  pressure, 
analyzed  mean  sea  level  pressure  and  terrain  height  values  were  used  with 
analyzed  geopotential  heights  to  interpolate  values  of  surface  pressure  for 
each  grid  point.  The  interpolation  was  based  on  the  assumption  that  geo¬ 
potential  height  is  a  quadratic  function  of  the  natural  logarithm  of  pres¬ 
sure,  following  the  practice  of  the  NMC.  The  cubic  spline  was  then  used  to 
obtain  values  of  q  at  even  multiples  of  50  mb  between  1000  mb  and  300  mb, 
at  psfc  £  1000  mb,  and  at  Psfc~  300.  In  cases  where  pgfc  >  1000,  qgfc  was 
set  equal  to  The  interpolation  and  integration  was  performed  in  the 

same  way  as  for  the  Level  II  data. 

To  perform  zonal  and  spherical  harmonic  analyses,  global  grid  point 
representations  of  the  observations  had  to  be  prepared.  In  choosing  a  meth¬ 
od  to  combine  analyses  and  observations  into  a  gridded  field,  it  was  impor¬ 
tant  to  minimize  the  modification  of  the  true  spectral  character  of  the  at¬ 
mospheric  field.  This  was  why  a  direct  interpolation  of  observations  onto 
grid  points  was  not  used,  since  most  polynomial-type  direct  interpolations 
have  a  bias  toward  certain  wave  numbers  (Yang  and  Shapiro,  197311) .  On  the 
other  hand,  the  Flattery  analysis  itself,  used  to  prepare  the  Level  III 
analyses,  also  includes  certain  biases  that,  without  attempted  correction, 
would  influence  the  harmonics  to  be  observed  here  in  the  results.  For  ex¬ 
ample,  the  Flattery  analysis  includes  a  truncation  at  zonal  wave  number  24, 
even  though  the  Level  III  analysis  values  are  available  at  2.5°  latitude- 
longitude  intervals,  which  would  allow  resolving  waves  as  small  as  zonal 
wave  number  72.  Furthermore,  the  Flattery  analysis  uses  12-hour  forecasts 
from  the  National  Meteorological  Center's  coarse-grid  hemispheric  model  as 
the  basis  for  the  first-guess  field.  The  Flattery  scheme  is  iterated  several 

11.  Yang,  C.  H.,  and  R.  Shapiro,  1973:  The  effects  of  the  observational 
system  and  the  method  of  interpolation  on  the  computation  of  spectra. 

J.  Atmos.  Sci..  4,  530-536. 


times,  beginning  with  the  expansion  in  terms  of  just  the  lowest  wave  numbers, 
and  including  higher  wave  numbers  with  successive  iterations.  (This  is  anal¬ 
ogous  to  the  successive  approximation  techniques  that  use  successively  larger 

scan  radii  in  the  selection  of  observations  to  influence  the  value  at  a  grid 

12 

point.)  In  addition,  chu  and  Parrish  (1977)  showed  that  a  simple  scheme 
that  transfers  data  to  grid  points  resulted  in  a  humidity  field  that  better 
represented  observations  than  the  Flattery  analysis  (described  by  Cooley, 
197413)  for  a  case  over  the  eastern  (J.  S.  While  admittedly  the  emphasis  in 
the  present  study  is  on  spectral  representation  of  the  observations  rather 
than  point-wise  accuracy,  the  nature  of  the  larger  point-wise  error  shown  in 
their  study  clearly  would  have  a  negative  impact  on  the  spectral  representa¬ 
tion.  Despite  these  shortcomings  of  the  Flattery  analysis,  it  has  the  favor¬ 
able  property  that  it  lacks  any  distinct  bias  toward  individual  groups  of 
waves,  since  it  is  based  on  an  expansion  of  the  data  in  unweighted  Hough 
basis  functions.  For  this  reason,  the  Flattery  analyses  were  used  as  the 
first  guess  field,  and  each  corresponding  set  of  observations  was  used  to 

correct  the  field  in  successive  scans  univariately  using  the  well  known  meth- 

14 

od  of  Cressman  (1959)  .  In  this  way,  corrected  fields  of  zonal  wind  (u) , 

meridional  wind  (v) ,  temperature,  geopotential  height  (Z) ,  specific  humidity, 
and  layer  precipitable  water  (W)  at  both  850  mb  and  500  mb  levels  were  pre¬ 
pared  for  each  of  the  14  observation  times  in  the  two  periods. 

To  measure  the  impact  of  the  corrections,  two  quantities  were  then  cal¬ 
culated  for  each  of  the  two  fields.  The  original  Flattery  analyzed  field 
(FG)  and  the  final  corrected  field  (CFG)  were  interpolated  bilinearly  to  the 
observation  locations  for  comparison  with  the  observed  values  (OB) .  Then 
the  following  two  root  mean  square  error  quantities  were  calculated  for  each 


field: 
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12.  Chu,  R.,  and  D.  Parrish,  1977:  Humidity  Analyses  for  Operational  Predic¬ 
tion  Models  at  the  National  Meteorological  Center.  Office  Note  140, 
National  Meteorological  Center,  NWS,  NQAA. 


Cooley,  D.  S.,  1974:  A  Description  of  the  Flattery  Global  Analysis 
Method  -  No.  1.  Tech.  Proc.  Bui.  No.  105,  NWS,  NOAA. 


14.  Cressman,  G.,  1959:  An  operational  objective  analysis  system.  Mon. 
Wea.  Rev.,  87,  367-374. 
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where  N  is  the  number  of  observations  used  in  each  field  for  correction  >f 
the  first  guess.  The  quantity  r  ■=  1  -  is  a  measure  of  the  im¬ 

pact  the  observations  had  in  modifying  the  analyzed  field  to  be  more  repre¬ 
sentative  of  the  observations;  the  closer  r  is  to  unity,  the  more  the  cor¬ 
rected  field  is  like  the  observed  field  with  respect  to  the  first  guess. 
Values  for  N  and  r  averaged  over  the  14  observation  periods  for  each  quan¬ 
tity,  level,  and  month  are  presented  in  Table  1. 

TABLE  1.  AVERAGE  NUMBER  OF  OBSERVATIONS  (N)  AND  ROOT  MEAN  SQUARE  ERROR 
RATIO  (r) 


DATE 

LEVEL 

ii 

v 

T 

Z 

a 

W 

850 

N 

810 

810 

608 

613 

587 

596 

JAN 

r 

.36 

.36 

.45 

.46 

.61 

.62 

500 

N 

718 

718 

624 

620 

545 

465 

r 

.32 

.33 

.38 

.45 

.60 

.62 

850 

N 

848 

848 

618 

616 

601 

629 

r 

.37 

.34 

.42 

.45 

.51 

.44 

JUL 

500 

N 

762 

762 

640 

634 

564 

554 

r 

.33 

.35 

.34 

.44 

.43 

.46 

In  most  cases,  the  table  shows  that  the  correction  had  a  somewhat  great¬ 
er  impact  on  T  and  Z  than  on  u  and  v,  with  a  still  greater  change  for  q  and 
W.  The  use  of  fewer  vertical  empirical  orthogonal  functions  in  the  Flattery 
humidity  analysis  might  be  suspected  of  creating  a  first  guess  humidity  that 
was  more  unlike  the  observations  than  the  first  guess  fields  for  the  other 
quantities.  However,  in  the  previous  study  (Gerlach,  19816)  it  was  shown 
that  the  first  three  vertical  modes,  the  ones  used  in  the  humidity  analysis, 
contain  80-90  percent  of  the  total  horizontal  variance  when  complete  radio¬ 
sonde  observations  were  analyzed  following  the  vertical  expansion  method  used 
in  the  Flattery  analysis.  Therefore,  it  is  unlikely  that  the  greater  change 
in  the  humidity  fields  can  be  explained  by  the  use  of  fewer  vertical  inodes  to 
create  the  first  guess  field.  Actually,  no  clear  explanation  can  be  given 
for  the  trends  in  the  r  values  in  Table  1.  However,  although  as  a  group 
they  show  that  substantial  point-wise  changes  result  from  tne  observational 
corrections,  the  corrected  fields  still  are  for  the  most  part  more  a  product 
of  the  first  guess  than  they  are  of  the  observations.  Apparently,  the  num¬ 
ber  of  observations  used  is  not  the  strongest  factor  in  determining  the 
amount  of  change  brought  about  by  the  corrections.  In  fact,  just  the  oppo¬ 
site  seems  to  be  the  case,  and  this  may  be  explained  by  the  fact  that  more 
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observations  available  in  the  Flattery  analysis  would  result  in  a  better 
representation  of  the  observations  by  the  first  guess.  Then  the  change 
brought  about  by  the  observations  would  be  expected  to  be  smaller.  The 
effect  on  these  results  from  the  interpolation  to  observation  sites  is  not 
considered.  A  better  method  would  have  been  to  use  the  coefficients  of  the 
Flattery  expansion  to  reproduce  Flattery  values  at  the  observation  sites, 
but  these  coefficients  were  not  available. 

III.  Harmonic  Analysis  of  Corrected  Fields 

Each  of  the  corrected  fields  was  analyzed  spectrally  using  both  zonal 
and  spherical  harmonic  expansions.  To  prepare  each  field  for  these  expan¬ 
sions,  the  departure  from  the  global  average  was  calculated  for  each  of  the 
grid  points  of  the  2.5°  latitude-longitude  global  grid.  The  global  average 
of  each  quantity  was  calculated  by  evaluating  the  integral 

in  . 

Ca]  =  h  J  A  ($)  cos  $  d$  (5) 

-tt/2 

using  a  Gaussian  quadrature  integration  in  the  form 

[»]  ■  H  £  iy.t  (6) 

where  <)>  represents  latitude,  A  represents  the  zonal  mean  of  the  quantity, 

(computed  as  the  arithmetic  average  of  the  144  points  around  each  latitude 

circle) ,  and  are  the  Gaussian  weights  corresponding  to  each  of  the  64 

—  o 

Gaussian  latitudes  used  in  the  study.  The  A  (4>)  values  at  every  2.5  lati¬ 
tude  were  interpolated  linearly  to  the  Gaussian  latitudes,  the  quadrature 
was  performed,  and  the  departure  value  A*  ($,X)  ■  A  ($,X)  -  [A]  was  calcu¬ 
lated  for  each  grid  point. 

Each  departure  field  was  expanded  along  each  2.5°  latitude  using  the 
harmonic  expansion  of  the  form 

M 

A*  ($,X)  ■  E,.  [a  ($)  cos  mX  +  b  (4>)  sin  mX]  (7) 

m*o  m  m 

where  X  is  longitude  and  m  is  the  zonal  wave  number  out  to  M  *  72.  The  coef¬ 
ficients  were  obtained  using  a  fast  Fourier  transform  algorithm. 

2  2  2 

The  zonal  mean  power  an  and  the  zonal  eddy  power  a  +  b  ,  m  =  1,  72  were 

o  ri  m 

calculated  for  each  latitude  for  each  quantity,  level,  and  time,  and  were 
then  averaged  over  time.  The  weekly  averaged  values  of  total  eddy  power  were 
normalized  in  the  form 
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(8) 


1“ 

m 


('(-) 


,  2  .2. 

{ a  i  h  I 
in  m 


£  {a2  +  b2} 
m=l  m  m 


where  the  brackets  {}  represent  the  weekly  averaged  quantity.  Thus  the  frac¬ 
tional  contribution  V  over  all  wave  numbers  for  a  particular  latitude  should 

m 

add  up  to  unity. 

Before  discussing  the  results,  a  comment  should  be  made  about  the  sensi¬ 
tivity  of  the  spectra  of  the  corrected  field  to  the  density  of  observations. 
With  the  addition  of  more  observations,  the  largest  changes  in  the  spectra 
would  be  seen  in  the  high  wave  numbers,  since  a  greater  density  of  observa¬ 
tions  would  allow  the  resolving  of  smaller  waves  in  the  fields  if  they  exist. 
The  long  waves  are  probably  already  well  represented  by  the  existing  less 
dense  data  coverage,  so  large  relative  changes  in  the  low  wave  numbers  would 
not  be  expected  by  the  addition  of  other  types  of  observations,  such  as  those 
from  satellites.  The  question  is  whether  the  present  data  may  produce  spec¬ 
tra  that  are  more  representative  of  the  actual  wave  structure  for  one  quan¬ 
tity  than  they  are  for  another.  In  this  regard,  for  example,  the  rawinsonde 
network  might  be  adequate  for  depiction  of  the  true  spectral  characteristics 
of  the  height  field  but  may  not  be  sufficiently  dense  to  resolve  important 
contributions  from  high  wave  numbers  in  the  q  field.  A  sensitivity  study  of 
the  effect  of  additional  observatons  on  the  various  spectra  was  not  included 
in  this  study,  but  it  would  be  useful  in  verifying  the  general  conclusions 
reached  in  the  present  project. 

2 

Values  of  weekly  averaged  zonal  mean  power  {a  }  and  zonal  eddy  power 
™  2  2  " 

£,  {a  +  b  }  are  shown  graphically  in  Pigs.  3a-d.  The  zonal  mean  power 
in  in 

curves  for  T,  Z,  q,  W  are  similar  in  that  they  all  show  the  maximum  positive 
departure  from  the  global  average  at  the  equator  and  maximum  negative  depar¬ 
ture  at  the  poles.  The  two  minima  in  the  curves  are  at  or  near  the  latitudes 
where  the  zonal  mean  of  the  quantity  is  equal  to  the  global  mean.  Zonal 
mean  power  curves  for  u  and  v  do  not  show  this  characteristic  variation. 

The  u  variation  reflects  the  latitudinal  variation  of  the  westerlies  at  tem¬ 
perate  latitudes  and  easterlies  in  tropical  latitudes.  The  zona lly-aver aged 
departure  from  the  global  mean  for  v  is  an  indicator  of  the  poleward-equator- 
ward  mass  flux  at  each  level,  a  small  quantity. 
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Fig.  3a.  Weekly  averaged  zonal  mean 
eddy  power  x — x — x  for  JAN  at  850  mb 


Pig.  3b.  Meekly  averaged  aortal  mean  power  - and  zonal 

eddy  power  x — x — x  for  JAN  at  500  mb. 
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Fig.  3d.  Weekly  averaged  zonal  mean  power  -  and  zonal 

eddy  power  x — x — x  for  JUL  at  500  mb. 
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2 

For  all  plots  of  {aQ}  for  T,  Z,  q,  and  W,  the  winter  hemisphere  minimum 
is  closer  to  the  equator  than  the  summer  hemisphere  minimum.  This  means 
that  the  zonal  average  of  each  quantity  decreases  to  the  global  average 
faster  with  latitude  in  the  winter  hemisphere  than  it  does  in  the  summer 
hemisphere.  This  makes  sense  since  the  average  winter  hemisphere  value  of 
T,  Z,  q,  and  W  is  less  than  the  summer  hemisphere  value.  When  the  locations 
of  the  minima  are  compared  for  a  given  month  at  a  given  level,  however,  it 
is  clear  that  both  the  summer  and  winter  hemisphere  minima  are  closer  to  the 
equator  for  q  and  W  than  they  are  for  T  and  Z.  When  the  zonal  averages  were 
plotted  as  a  function  of  latitude  for  T,  Z,  and  q  for  1/15/78  at  00Z,  q 
showed  a  much  more  rapid  decrease  with  latitude  in  both  hemispheres  than  did 
the  other  two  quantities.  This  steeper  meridional  gradient  accounts  for  the 
closer  position  of  the  zonally  averaged  global  mean  departures  to  the  equator 
for  q  and  W.  The  latitudinal  difference  between  the  summer  and  winter  hemis¬ 
phere  minima  decreases  significantly  with  height  for  q  and  W,  mostly  due  to 
the  movement  of  the  winter  hemisphere  minimum  to  lower  latitudes.  This 
shift  is  not  nearly  as  pronounced  in  T  and  Z.  This  indicates  that  the  mag¬ 
nitude  of  the  winter  hemisphere  meridional  gradient  of  the  moisture  variables 
must  increase  with  height,  unlike  that  of  the  temperature  and  height  fields. 

The  total  zonal  eddy  power  curves  for  T  and  Z  show  a  bimodal  character, 
with  maxima  at  mid-latitudes  and  minima  at  the  equator  and  both  poles.  This 
is  also  true  to  a  lesser  extent  of  the  u  and  v  distributions.  This  is  not 
surprising  due  to  the  large  amplitude  eddies  which  form  in  and  propagate  in 
the  westerlies  at  mid- latitudes.  In  contrast  to  this  is  the  latitudinal 
distribution  of  zonal  eddy  power  for  q  and  W,  which  shows  a  maximum  at  the 
equator  and  minima  at  the  poles.  Apparently,  the  amplitudes  of  the  depar¬ 
tures  from  the  zonal  mean  over  all  length  scales  are  greater  in  the  tropics 
than  those  in  mid-latitudes.  Two  possible  explanations  for  this  are  (a) 
the  restriction  of  the  amplitude  of  q  (and  thus  W)  about  the  zonal  mean  re¬ 
sulting  from  the  saturation  specific  humidity  qg  which  is  lower  in  the  mid¬ 
latitudes  than  it  is  in  the  tropics,  and  (b)  the  larger  magnitude  of  sta¬ 
tionary  departures  from  the  zonal  mean  of  q  and  W  (due  to  land-ocean  distri¬ 
bution)  in  the  tropics.  Dealing  with  the  latter  possibility  first,  one 
could  measure  the  contribution  of  this  effect  by  removing  the  time-averaged 
component  from  each  of  the  fields  and  performing  the  same  zonal  harmonic 
analysis  of  the  time  departure  fields,  then  comparing  the  resulting  total 
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eddy  power  curves  with  those  of  Pig.  3.  However,  the  one-week  time  period 
used  in  each  case  would  probably  not  be  a  sufficient  length  of  time  to  as¬ 
sure  that  the  time-average  is  representative  of  climatological  conditions. 
Using  the  climatological  values  of  W  in  Fig.  2  along  with  vertically  inte¬ 
grated  values  of  g  from  Oort  and  Rasmusson  (1971)^  for  a  yearly  average, 

o 

the  largest  apparent  departure  from  the  zonal  mean  south  of  20  N  was  about 

_2  M  2  0 

2.0  x  10  kg  m  ,  compared  to  a  departure  of  1.2  x  10  kg  m  at  40  N.  In  fact, 

the  largest  departure  of  the  time-averaged  field  from  the  zonal  mean  occurs 

at  30°N,  most  likely  resulting  from  the  Sahara.  This  large  departure  shows 

up  at  30°N  on  the  zonal  eddy  power  plots  in  Figs.  3a  and  3c  for  q  and  W, 

expecially  in  the  JUL  case.  The  fact  that  at  least  qualitatively  the  time 

averaged  departures  below  20°N  are  not  significantly  larger  than  mid-latitude 

zonal  departures  and  that  the  30°N  maximum  shows  up  as  just  a  prominence  on 

already  large  low  latitude  values  of  eddy  power  for  q  and  W  suggests  that 

factor  (b)  is  not  primarily  responsible  for  the  general  shape  of  the  q  and 

W  curves. 

To  consider  the  case  for  the  first  explanation,  consider  another  effect 
of  a  lower  qg  value  resulting  from  lower  temperatures.  By  comparing  the 
right  hand  sides  of  the  zonal  power  curves  for  JAN  with  the  corresponding 
section  of  curve  for  JUL  (the  Northern  Hemisphere  sections  are  compared  due 
to  their  higher  observation  density),  it  is  clear  that  for  u,  v  and  espe¬ 
cially  T,  Z,  the  total  zonal  eddy  power  (area  under  the  curve  for  NH)  is 
greater  for  the  winter  hemisphere  than  it  is  for  the  summer  hemisphere. 

Again,  this  is  not  surprising  because  the  amplitudes  of  large-scale  wave  dis¬ 
turbances  are  clearly  larger  in  the  winter  hemisphere  due  to  greater  baro- 
clinicity.  However,  the  q  and  W  power  curves  show  greater  zonal  eddy  power 
in  the  summer  hemisphere.  It  appears  that  higher  values  for  the  q  and  W 
curves  are  favored  in  warm  regions  and  lower  values  in  the  cooler  regions. 
Thus,  it  could  be  that  the  lower  temperatures  tend  to  reduce  the  total  var¬ 
iance  about  the  zonal  mean  because  of  the  limiting  value  of  qg  which  de¬ 
creases  with  decreasing  temperature.  Clear  support  for  this  argument  can  be 
gained  from  Table  2.  The  corrected  fields  of  T  for  00GMT  on  1/15/78  for 
both  levels  were  used  to  generate  fields  of  saturation  specific  humidity  qg. 
Then  the  q  fields  for  the  same  time  and  date,  whose  zonal  power  curves  look 
very  much  like  the  weekly  average  curves  in  Fig.  3,  were  used  with  the  qg 
fields  to  form  RH  =  q/qg  at  each  grid  point.  The  Northern  Hemisphere  aver¬ 
age  values  show  that  while  temperature  increases  by  only  4  percent  between 
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winter  and  summer,  the  saturation  specific  humidity  q  ,  acting  as  an  upper 

s 

limit  for  q,  increases  by  73  percent  at  850  mb  and  59  percent  at  500  mb. 

For  850  mb,  this  translates  to  a  maximum  allowable  departure  from  the  aver¬ 
age  in  q  of  2.6  g/kg  in  the  winter  compared  to  4.9  g/kg  in  the  summer.  The 
hemispheric  average  relative  humidity  is  more  like  the  temperature,  in  that 
there  is  very  little  change  between  the  seasons.  This  suggests  that  the 
limiting  magnitude  of  the  departure  about  the  mean  for  relative  humidity  is 
not  nearly  as  sensitive  to  temperature  as  it  is  for  specific  humidity. 


TABLE 

2 .  NORTHERN 

HEMISPHERE 

AVERAGES 

q(g  kg”1) 

T(°K) 

qs(g  kg”1) 

RH 

850 

1/15/78 

5.35952 

276.49600 

8.00103 

.66985 

850 

7/20/78 

8.95516 

288.18456 

13.83457 

.64730 

500 

1/15/78 

1.23598 

254.71750 

2.67001 

.46291 

500 

7/20/78 

2.13636 

267.28384 

4.25537 

.50204 

Further  evidence  for  the  importance  of  the  greater  limiting  influence 

of  qg  can  be  gained  from  Figs.  4a-d,  the  zonal  eddy  power  curves  of  RH  for 

1/15/78  and  7/20/78.  Two  facts  can  be  noticed:  (a)  the  total  power  in  the 

winter  hemisphere  is  as  large  as  that  for  the  summer  hemisphere,  and  (b)  the 

equatorial  maximum  seen  in  the  q  and  W  curves  is  not  evident  in  the  RH 

curves.  In  fact,  there  is  indication  of  a  bimodal  character  in  the  curves 

as  in  the  T  and  Z  curves,  especially  at  850  mb.  The  table  values  and  curves 

would  appear  to  explain  that  the  variance  of  q  in  colder  regions  is  held 

down  by  the  lower  value  of  q  ,  and  that  relative  humidity  is  not  so  strongly 

s 

affected  by  this  temperature  limitation. 

Values  of  F  ,  the  normalized  (fractional)  zonal  eddy  power,  are  contour- 
m 

ed  in  1°910  form  (solid  contours)  in  Figs.  5a-d.  Only  the  values  for  the 
first  32  zonal  wave  numbers  were  included  in  the  plots,  since  power  contri¬ 
butions  from  higher  wave  numbers  were  negligibly  small.  The  contours  are 
labeled  by  the  power  of  10  corresponding  to  the  fractional  contribution  to 
the  total  eddy  power  at  a  particular  latitude  by  a  particular  wave  number. 

For  example,  a  point  lying  to  the  left  of  the  -2  contour  indicates  that  that 
wave  number  contributes  at  least  one  percent  of  the  total  power  at  the  cor¬ 
responding  latitude.  The  daBhed  curves  represent  equatorial  wave  number  — 
any  point  lying  along  a  given  curve  represents  a  constant  wavelength  equal 
to  the  wavelength  corresponding  to  that  number  of  waves  around  the  equator. 
The  dashed  curves  were  included  because  it  was  evident  that  the  normalized 
power  contours  generally  followed  wavelength  rather  than  local  wave  number. 
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total  power 


Fig.  4a.  Total  zonal  power  for  the  departures  from  the  zonal 
average  for  the  ratio  q/qg.  Both  q  and  T  used  in  the  calcu¬ 
lations  are  the  corrected  field  values. 
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Fig.  4b.  Total  zonal  power  for  tha  departures  from  the  zonal  average 
for  the  ratio  q/q  .  Both  q  and  T  used  in  the  calculations  are  the 
corrected  field  vllues. 
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Fig.  4c.  Total  tonal  power  for  the  departure*  front  the  eonal  average 
for  the  ratio  q/q  .  Both  q  and  T  used  in  the  calculation*  are  the 
corrected  field  vllua*. 
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RH  500  7/20/78 


Fig.  4d.  Total  zonal  power  for  the  departures  from  the  zonal  average 
for  the  ratio  q/q  .  Both  q  and  T  used  in  the  calculations  are  the 
corrected  field  values. 
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Fig.  5a.  Meekly  averaged  normalized  zonal  eddy  power  for  JAN 
at  850  ab.  Contour  labelB  are  powers  of  10  corresponding  to 
the  fractional  contribution  of  each  wave  nunber.  Dashed  curves 
represent  lines  of  constant  equatorial  wave  number,  given  by 
the  contour  labels. 


Pig.  5c.  Weekly  averaged  normalized  zonal  eddy  power  for  JUL 
at  850  mb . 
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This  is  especially  true  of  the  -2  contour.  Notice  the  truncation  of  frac¬ 
tional  power  at  wave  number  24  as  evidenced  by  the  vertical  sections  of  con¬ 
tour  along  this  wave  number.  This  is  a  remnant  from  the  truncation  of  zonal 
wave  number  24  imposed  in  the  Flattery  analysis.  In  comparing  the  g  and  W 
curves  with  those  for  the  other  four  quantities,  one  observes  that  for  v,  g, 
and  W,  the  power  is  more  spread  out  toward  higher  wave  numbers  than  is  the 
case  with  the  u,  T  and  Z  plots.  Again,  the  -2  contour  shows  this  most  clear¬ 
ly.  For  example,  the  -2  (1%)  contour  lies  generally  along  the  eguatorial 
wave  number  20  curve  for  g,  whereas  for  T  it  follows  the  wave  number  16  curve, 
especially  in  the  Northern  Hemisphere.  The  T  and  Z  curves  tend  to  show  a 
restriction  of  power  to  lower  wave  numbers  in  mid-latitudes,  and  more  spec¬ 
tral  spreading  in  the  tropics.  Such  a  pattern  is  not  evident  in  the  g  and 
W  plots,  in  which  the  1%  contour  generally  follows  the  dashed  constant  wave¬ 
length  curves  at  all  latitudes. 

Table  3  gives  a  guantitative  comparison  of  the  zonal  eddy  power  distri- 

o 

bution  by  wave  number  bands  at  40  N.  At  this  latitude,  the  shift  of  zonal 
eddy  power  to  higher  wave  numbers  with  height  for  g  and  W  is  clearly  evi¬ 
dent  for  JAN,  but  is  totally  absent  in  the  JUL  case.  Just  the  opposite  is 

true  for  T  and  Z,  in  which  no  power  shift  with  height  is  evident  in  the  JAN 

case  but  is  guite  pronounced  in  the  JUL  case.  This  same  shift  is  evMent 

to  a  lesser  extent  for  u.  The  power  distribution  for  v  remains  consistent 
for  both  cases  and  levels,  with  the  vast  majority  of  the  zonal  eddy  power 
lyir  .n  the  middle  band. 


TABLE  3.  FRACTIONAL  CONTRIBUTIONS  TO  TOTAL  EDDY  POWER  IN  ZONAL  HARMONICS 
SPECTRUM  FOR  40°N  FROM  THREE  WAVE  NUMBER  BANDS 


DATE 


JAN 


JUL 


LEVEL 

WAVE  NUMBER 

u 

v 

T 

Z 

q 

W 

BAND 

■ 

_ 

_ 

_ 

1-2 

.46 

.08 

.55 

.61 

.46 

.50 

850 

3-10 

.47 

.72 

-39 

.37 

.41 

.39 

11-20 

06 

.17 

.05 

.02 

.11 

.09 

1-2 

.42 

.04 

59 

.57 

.31 

.31 

500 

3-10 

.52 

.77 

.35 

.41 

.52 

.53 

11-20 

.06 

.16 

.05 

.02 

.14 

.13 

1-2 

.36 

.13 

.63 

.71 

.25 

.31 

850 

3-10 

.50 

.65 

.31 

.28 

.53 

.51 

11-20 

.11 

.18 

.05 

.01 

.17 

.13 

1-2 

.16 

.04 

.32 

.35 

.24 

.31 

500 

3-10 

.67 

.69 

.55 

.59 

.53 

.50 

11-20 

.15 

.24 

.11 

.06 

.19 

.15 
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After  the  zonal  harmonic  analysis,  each  global  departure  field  was  ex¬ 


panded  in  a  spherical  harmonic  expansion  of  the  form 

M  m+N 

A  *(<p,\)  =  I  Z 

m=0  n=m 


r  n\  ,  _  m  ,  .  -«  m  . 

La  cos  mX  +  B  sin  mX  ]  P  (sm  *) 
n  n  n 


(9) 


where  N  represents  the  truncation  limit  of  the  quantity  n-m,  which  for  the 
associated  Legendre  function  P™  (sin  4>)  represents  the  number  of  times  the 
function  has  a  value  of  zero  between  the  poles.  Thus  n-m  can  be  thought  of 
as  a  type  of  meridional  wave  number  as  a  counterpart  to  the  zonal  wave  num¬ 
ber  m,  whose  truncation  limit  is  set  at  M  =  32  for  this  expansion.  To  eval¬ 
uate  the  coefficients  am  and  Bm,  the  A*  values  at  the  grid  points  were  in- 

n  n 

terpolated  linearly  to  the  Gaussian  latitudes,  and  a  fast  Fourier  transform 

was  used  to  evaluate  the  coefficients  y  and  6  of  the  zonal  expansion 

mk  mk 


M 

A*  (<t>,  ,X)  =  Z  Y  cos  mX  +  6  sin  mA 

k  „  m.  m, 

m*0  k  k 


(10) 


where  k  represents  the  indices  of  the  Gaussian  latitudes  along  which  this 


out. 

and  Y 
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n=m 

m+N 

5 

=  l 

8m 

Pm 

n 

n 

n=m 

(11) 


where  fi,  =  sin  <)>,  .  Using  the  orthonormality  of  Pm  with  values  of  P®  that 
K  k  nr. 

normalize  to  unity,  we  have 


m 

r 

m 

i  * 

1  y 

P 

n 

J ,  m. 

n 

-1  k 

1 

c  , 

m 

f  6 

P 

n 

J  ■.  m , 

n 

-1  k 

(12) 


These  integrals  were  evaluated  using  the  Gaussian  quadrature  technique. 

For  M  »  32  and  64  Gaussian  latitudes,  the  highest  order  Pm  that  can  be  in- 

32  n 

tegrated  exactly  is  P  ,  so  N  »  31  was  chosen  to  complete  the  rhomboidal 

o3 

truncation  of  the  spherical  harmonic  expansion.  Once  these  numerical  in¬ 
tegrations  were  carried  out,  the  weekly  averaged  spherical  power  was  used 
to  calculate  the  normalized  spherical  harmonic  power: 
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(13) 


Figs  6a-d  display  A™  plotted  in  log1Q  ^orm  on  an  n"m  vs*  m  P^ot>  The 
contour  labels  are  the  powers  of  10  corresponding  to  the  fractional  contri¬ 
butions  at  each  (n-m,m)  ordered  pair  of  modes.  The  sum  of  the  fractional 
contribution  from  each  ordered  pair  over  the  entire  plot  is  unity.  Here, 
the  -3  and  -5  "ontours  alone  are  shown  for  clarity  —  they  correspond  to  a 
0.1  and  0.001  percent  contribution  to  the  total  spherical  power,  respective¬ 
ly.  The  24  wave  number  truncation  is  somewhat  evident  in  the  plots  as  was 
the  case  with  the  zonal  power  plots.  In  these  plots,  the  largest  spectral 
spreading  seems  to  be  occurring  in  the  wind  components.  Both  contours  on 
the  v  diagrams  lie  at  relatively  high  wave  numbers  in  both  directions,  while 
the  u  curves  are  spread  out  in  primarily  the  meridional  direction.  This  is 
certainly  consistent  with  the  plot  of  the  zonally  averaged  departure  of  u 
in  Fig.  3  which  showed  a  predominance  for  several  waves  from  pole  to  pole. 
The  wind  component  plots  show  very  little  variation  with  height  in  compari¬ 
son  with  the  other  four  quantities.  While  the  contours  on  the  q  and  W  plots 
lie  at  higher  wave  numbers  in  both  directions  at  500  mb  than  at  850  mb  for 
both  months,  just  the  opposite  is  true  for  Z  and  to  a  lesser  degree,  T. 

Also,  for  any  given  month  and  height,  it  is  clear  that  T  and  Z  have  a 
greater  percent  of  their  power  in  lower  wave  numbers  than  do  the  other  four 
quantities.  In  general,  it  appears  that  T  and  Z  are  most  restrictive  in 
concentrating  power  in  low  wave  numbers,  q  and  W  are  less  so,  and  u  and  v 
hold  a  greater  share  of  their  spectral  power  in  higher  wave  numbers. 

These  findings  are  illustrated  in  another  way  in  Figs.  7a-b.  Here  cum¬ 
ulative  fraction  of  the  power  contributed  by  all  of  the  components  that  have 
both  n-m  and  m  less  than  or  equal  to  the  value  of  "wave  number"  on  the  ab¬ 
scissa  is  plotted  against  wave  number.  Thus,  for  example,  the  value  plotted 
on  the  figure  at  wave  number  2  represents  the  sum  of  the  fractional  power 
contributions  from  all  wave  number  combinations  such  that  both  m  and  n-m  are 
less  than  or  equal  to  2.  The  large  jumps  on  the  graphs  at  "wave  number"  =  2 
show  the  dominance  of  the  n-m  ■*  2,  m  •  0  component  in  the  T,  Z,  q,  and  W 
harmonics.  While  this  harmonic  is  not  nearly  as  dominant  at  850  mb  as  it 
is  at  500  mb  for  Z,  the  dominance  is  nearly  the  same  at  both  heights  for  T, 
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Fig.  6a.  Weekly  averaged  normalized  spherical  eddy  power 
for  JAN  at  850  mb.  Contour  labels  are  powers  of  10  corres¬ 
ponding  to  the  fractional  contribution  of  each  wave  number 
(n-m,m)  ordered  pair.  The  value  n— m  is  the  number  of  zeros 
the  curve  passes  through  between  the  poles. 
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*W=^lywaVera9ed  normalized  spherical  eddy  power 
for  JAN  at  500  mb. 


I 


I 
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Fig.  6d.  Weekly  < 
for  JU1.  «t  500  mb 
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Fig.  7a.  Cumulative  normalized  eddy  power  for  JAM  at  850  mb 

-  and  500  mb  x — x — x .  "Wave  number"  on  the  abscissa  refer 

to  equal  valuee  of  n-m  and  m  on  Figs.  6a-d  over  which  the  f  ra'.' 
tional  contributions  were  summed  to  produce  these  plots. 


and  the  q  and  W  plots  reveal  a  greater  dominance  in  this  harmonic  at  850  mb. 
Primarily  because  of  this  and  the  fractional  contributions  up  to  about  "wave 
number”  12,  a  rhomboidal  truncation  set  at  12  in  an  analysis  would  clearly 
include  essentially  all  of  the  power  for  T  and  Z,  whereas  a  significant  por¬ 
tion  would  be  missing  for  q  and  W,  especially  at  500  mb.  Such  a  truncation 
would  create  an  even  more  serious  distortion  in  an  attempt  to  represent  the 
wind  field,  since  a  good  share  of  its  power  contribution  lies  in  wave  num¬ 
bers  greater  than  10.  This  is  revealed  more  quantitatively  in  Table  4,  in 
which  the  fractional  spherical  hormonics  contributions  from  three  "wave 
number"  bands  are  listed. 

TABLE  4.  CONTRIBUTIONS  TO  TOTAL  EDDY  POWER  IN  SPHERICAL  HARMONICS  SPECTRUM 


FROM  THREE 

WAVE  NUMBER 

BANDS 

DATE 

LEVEL 

WAVE  NUMBER 

u 

V 

T 

Z 

q 

W 

BAND 

1-2 

.13 

.03 

.82 

.60 

.65 

.66 

850 

3-10 

.70 

.65 

.14 

.37 

.25 

.25 

11-20 

.15 

.27 

.03 

.03 

.08 

.07 

JAN 

1-2 

.07 

.02 

.84 

.81 

.40 

.44 

500 

3-10 

.77 

.72 

.14 

.17 

.41 

.40 

11-20 

.14 

.22 

.02 

.01 

.16 

.14 

1-2 

.11 

.04 

.81 

.54 

.63 

.62 

850 

3-10 

.67 

.69 

.15 

.41 

.24 

.25 

11-20 

.19 

.23 

.04 

.04 

.10 

.09 

JUL 

1-2 

.11 

.01 

.85 

.81 

.37 

.41 

500 

3-10 

.70 

.75 

.12 

.17 

.41 

.39 

11-20 

.17 

.21 

.02 

.02 

.19 

.16 

A  significant  fraction  of  the  power  lies  in  the  11-20  band  for  u  and  v.  The 
fraction  from  this  band  is  less  for  the  moisture  quantities,  and  almost  in¬ 
significant  for  T  and  Z.  The  table  also  illustrates  the  relative  constancy 
of  u,  v,  and  T  between  season  and  height,  and  the  relative  shift  of  power 
with  height  of  Z  compared  to  q  and  W.  Finally,  it  shows  that  wave  number 
1-2  contributions  for  Z  are  slightly  less  than  for  q  and  W  at  850  mb  but 
rise  to  be  equal  to  T  at  500  mb. 

The  dominant  harmonic  for  T,  Z,  q,  and  W  is  the  n-m  =  2,  m  =  0  harmonic 
which  is  the  zonal  mean  of  each  quantity  which  has  a  maximum  at  the  equator 
and  minima  at  the  poles.  The  fractional  power  contribution  from  this  har¬ 
monic  is  given  in  xable  5,  along  with  the  contribution  from  the  n-m  *=  1, 
m  *  0  harmonic.  The  results  show  that  while  the  n-m  ■  2  harmonic  is  clearly 
the  dominant  one  in  all  cases,  the  n-m  *  1  (maximum  at  one  pole,  minimum  at 
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the  other)  becomes  larger  at  the  expense  of  the  n-m  =  2  component  in  JUL  in 
T  and  Z.  For  q  and  W,  the  seasonal  decrease  of  the  n-m  =  2  contribution  is 
not  as  large,  and  it  is  not  the  n-m  *  1  component  that  appreciably  benefits 
from  its  decrease.  The  probable  explanation  for  the  greater  n-m  *  1  contri¬ 
bution  in  the  summer  from  T  and  Z  is  the  higher  hemispheric  summer  tempera¬ 
ture  in  one  hemisphere  over  the  other.  When  zonally  averaged  temperatures 
for  1/15/78  and  7/20/78  (both  at  OOGMT,  850  mb)  were  plotted  as  a  function 
of  latitude,  (see  Figs.  8a, b)  the  1/15/78  case  depicted  much  greater  sym¬ 
metry  about  the  equator  than  the  7/20/78  plot.  The  7/20/78  case  clearly 
shows  a  combination  of  both  a  wave  number  1  and  2  component  in  it  (Northern 
Hemisphere  latitude  values  substantially  greater  than  Southern  Hemisphere 
counterparts,  and  absolute  maximum  shifted  northward)  whereas  in  the  1/15/78 
plot  there  is  no  clear  dominance  of  the  Southern  Hemisphere  values  over  cor¬ 
responding  Northern  Hemisphere  values.  The  same  was  true  at  500  mb  for  both 
T  and  Z  (not  shown) .  As  a  direct  result  of  this,  there  is  somewhat  less  sym¬ 
metry  about  the  equator  in  the  zonally  averaged  q  values  for  7/20/78  than 
there  is  for  1/15/78  (see  Figs.  8c, d) .  Since  these  comparisons  are  based  on 
only  one  observation  time  in  each  season,  the  results  are  not  conclusive. 
However,  on  the  basis  of  these  plots  it  appears  that  the  summer  zonally  aver¬ 
aged  temperature  is  higher  in  the  Northern  Hemisphere  than  in  the  Southern 
Hemisphere.  A  larger  land  to  ocean  ratio  in  the  Northern  Hemisphere  is  the 
probable  explanation  for  this,  due  to  greater  atmospheric  heating  over  land. 


TABLE  5. 


FRACTIONAL  CONTRIBUTION  FROM  THE  ZONAL  AVERAGE,  MERIDIONAL  WAVE 
NUMBER  1  AND  2  HARMONICS 


DATE 

JAN 

JUL 


LEVEL 

n-m 

T 

Z 

a 

W 

850 

1 

.08 

.01 

.05 

.05 

2 

.67 

.50 

.54 

.54 

500 

1 

.07 

.01 

.04 

.04 

2 

.73 

.76 

.31 

.36 

850 

1 

.25 

.13 

.09 

.08 

2 

.50 

.26 

.44 

.44 

500 

1 

.26 

.27 

.03 

.03 

2 

.57 

.53 

.25 

.28 

The  values  for  weekly  averaged  total  power,  the  denominator  in  the  ex- 

IB 

pression  for  A  ,  are  given  in  Table  6.  Because  the  values  are  dimensional 
n 

and  because  of  the  normal  variation  with  height  expected  in  each  variable, 
a  comparison  of  variation  in  season  only  can  be  made  between  the  variables. 


48 


TUN 


LA1ITU0C 


Fig.  8b.  Zonal ly  averaged  temperature  for  7/20/78  00GMT  at  850  mb 
computed  from  corrected  field. 


Fig.  8c.  Zonally  averaged  specifi' 
computed  for  corrected  field. 


LHTITU1C 


Fig.  8d.  Zonal ly  averaged  specific  tumidity  for  7/20/78  00GMT  at  850  mb 
computed  from  corrected  field. 


With  the  exception  of  v  at  850  mb,  all  of  the  quantities  except  q  and  W  show 
a  decrease  in  total  power  from  JAN  to  JUL.  The  zonal  eddy  power  was  greater 
in  the  winter  hemisphere,  and  it  is  likely  that  the  larger  number  of  obser¬ 
vations  in  the  Northern  Hemisphere  resulted  in  a  somewhat  more  irregular  cor¬ 
rected  field  in  data  rich  areas  compared  to  the  smoother  corrected  field  in 
data  sparse  areas,  contributing  to  the  larger  power.  The  lower  zonal  eddy 
power  for  the  moisture  variables  in  the  winter  hemisphere  (see  Fig.  3)  is 
probably  offset  by  the  greater  variation  (due  to  greater  data  density)  in 
ti.e  Northern  Hemisphere  moisture  fields  compared  to  that  of  the  Southern 
Hemisphere,  resulting  in  very  little  variation  in  total  power  with  season. 


TABLE  6. 

TOTAL  POWER  FOR  SPHERICAL  HARMONICS 

u  (m/s)2 

v  (m/s)2 

T  (°K) 2 

Z  (gpm) 2 

q  (g/kg)2 

W  (cm) 2 

JAN 

850 

200.7 

100.6 

349.0 

38103.4 

45.1 

4.6 

JAN 

500 

498.1 

250.5 

297.1 

179664.2 

6.7 

1.0 

JUL 

850 

167.8 

103.5 

251.9 

23104.6 

43.0 

4.1 

JUL 

500 

353.4 

202.1 

226.9 

120152.9 

5.8 

0.9 

IV.  Impact  of  Results  on  Objective  Analysis  Methods  for  Moisture 

It  is  evident  that  the  observed  moisture  field  has  certain  spectral  char¬ 
acteristics  that  make  it  uniquely  different  from  the  motion  and  mass  fields. 

In  the  preparation  of  initial  conditions  for  a  numerical  weather  prediction 
model,  any  gridded  representation  of  the  moisture  field  should  possess  these 
same  spectral  characteristics  in  order  to  provide  an  accurate  initial  condi¬ 
tion  for  a  global  model.  In  fact,  these  characteristics  can  also  be  used  to 
assess  the  accuracy  of  the  model  output  results  as  well.  As  mentioned  in 
the  beginning  of  this  report,  the  importance  of  analysis  accuracy  in  the 
model  forecast  will  depend  on  the  model  itself.  However,  the  more  faithful 
the  model  is  to  the  simulation  of  natural  processes,  the  more  it  will  de¬ 
pend  on  realistic  initial  conditions  to  generate  good  forecasts,  especially 
for  short-term  forecasts.  With  all  of  the  present  emphasis  on  forecast  mod¬ 
el  improvement,  it  is  certainly  timely  to  focus  attention  on  accuracy  in  the 
production  of  model  initial  conditions. 

Within  the  spectral  range  that  can  be  accurately  resolved  by  the  rawin- 
sonde  network,  the  spectral  results  of  this  study  are  an  accurate  representa¬ 
tion  of  the  spectral  characteristics  of  the  global  moisture  distribution. 

They  therefore  can  be  used  as  a  standard,  especially  in  the  lower  (less  than, 
say,  wave  number  20)  waves,  against  which  the  representations  from  any  analy- 
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sis  scheme  should  be  tested  for  their  accuracy.  The  zonal  and  spherical 
harmonic  analyses  of  the  gridded  fields  produced  by  the  analysis  scheme 
should  demonstrate  the  same  relative  features  as  those  that  are  shown  in 
this  study.  Any  major  differences  would  point  to  particular  problems  that 
may  exist  in  the  analysis  scheme  and  its  ability  to  represent  the  fields 
spectrally. 

The  faithfulness  of  analysis  schemes  in  the  spectral  representation  of 
the  observed  field  may  or  may  not  be  influenced  by  the  "closeness  of  fit" 
of  the  gridded  values  to  the  observations.  For  example,  if  the  amplitudes 
of  the  spectral  components  were  in  good  agreement  with  those  of  the  observa¬ 
tions  but  the  two  fields  were  out  of  phase,  the  analysis  would  test  out  well 
spectrally  but  would  be  a  failure  in  physical  space.  Both  a  "closeness  of 
fit"  in  physical  space  and  a  harmonic  analysis  in  spectral  space  must  be  per¬ 
formed  to  determine  if  the  analyzed  field  is  a  suitable  representation  of  the 
observations.  Both  tests  were  conducted  in  this  study  using  the  Flattery  an¬ 
alysis  as  the  test  analysis  scheme.  To  test  for  closeness  of  fit  in  physical 
space,  the  quantity 


1 

N 


2 

Z  (F*  -  OB*) 

.  ,  1  l 

1*1 


2 

Z  (OB*) 

i-1 


>5 


(14) 


was  calculated  using  the  Flattery  values  (F)  and  observations  (OB)  for 
1/15/78  OOGMT  and  7/16/78  OOGMT.  The  Flattery  analysis  values  were  inter¬ 
polated  bilinear ly  to  the  observation  sites,  their  global  mean  was  removed, 
and  they  were  compared  with  the  departures  of  the  N  observed  values  (de¬ 
partures  from  the  Flattery  global  mean,  found  to  be  very  close  to  global 
mean  calculated  from  corrected  fields)  for  u,  v,  T,  Z,  and  q.  The  smaller 
the  value  of  this  quantity,  the  better  the  overall  fit  of  the  analysis  field 
to  the  observations.  Table  7  lists  the  results  of  these  calculations.  The 
values  for  T  and  Z  are  generally  the  best,  except  for  Z,  1/15/78,  at  850  mb. 
The  accuracy  of  the  analysis  of  u  and  v  seems  to  improve  with  height,  but 
in  all  cases  it  is  not  as  good  as  the  accuracy  of  T  and  Z.  The  accuracy  of 
the  analysis  of  q  seems  to  be  more  seasonal  than  the  rest  of  the  quantities, 
with  accuracy  decreasing  from  winter  to  summer  by  about  the  same  amount  at 
both  levels.  It  appears  on  the  basis  of  these  tests  that  q  is  a  little 
more  accurately  analyzed  than  are  u  and  v,  but  less  accurately  than  T  and  Z. 
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TABLE  7.  NORMALIZED  RMS  ERROR  OF  FLATTERY  ANALYSIS 


1/15/78 

7/16/78 

850 

u 

.51 

.59 

V 

.54 

.76 

T 

.18 

.23 

Z 

.36 

.16 

q 

.30 

.55 

500 

u 

.31 

.37 

V 

.34 

.33 

T 

.13 

.20 

Z 

.09 

.13 

q 

.30 

.58 

A  direct  spectral  comparison  of  the  Flattery  analyses  with  the  observa¬ 
tions  is  not  possible  due  to  the  spatial  irregularity  of  the  observations. 

The  best  that  can  be  done  is  a  comparison  of  the  Flattery  power  spectra 
with  those  of  the  corrected  fields,  which  was  suggested  earlier  as  a  basis 
for  spectral  accuracy  of  the  analysis.  Such  a  comparison  is  given  for  the 
January  weekly  averages  in  Figs.  9a, b.  The  solid  curves  are  the  base  10 
logarithm  of  the  fractional  spherical  harmonic  power  for  the  Flattery  fields, 
while  the  dashed  curves  are  a  reproduction  of  the  curves  for  the  corrected 
fields  from  Fig.  6.  Except  for  larger  contributions  to  the  total  power  from 
higher  wave  numbers,  the  figures  do  not  show  very  much  difference  in  spec¬ 
tral  configuration  between  the  Flattery  and  corrected  fields.  In  fact,  ex¬ 
cept  for  the  case  of  Z,  and  to  a  lesser  extent  T,  it  appears  that  the  change 
brought  about  by  the  corrections  consisted  primarily  in  adding  back  in  the 
high  wave  number  variation  that  was  removed  in  the  wave  number  24  truncation 
of  the  Flattery  analysis.  This  suggests  that,  spectrally,  the  major  differ¬ 
ence  between  the  two  fields  is  in  the  amplitudes  of  the  shortest  waves.  The 
total  power  was  not  changed  significantly  by  the  corrections,  either.  For 
the  weekly  averaged  January  values,  the  total  power  increased  from  1  percent 
for  T  at  500  mb  to  12  percent  for  v  at  850  mb.  Thus,  the  apportionment  of 
power  among  the  wave  numbers  and  the  total  spherical  power  are  very  close 
to  the  same  in  both  fields. 

In  this  case,  there  is  favorable  spectral  agreement  between  the  cor¬ 
rected  fields  and  Flattery  analysis  fields,  but  according  to  Table  7  there 
are  appreciable  errors  when  the  Flattery  fields  are  compared  to  the  obser¬ 
vations.  This  apparent  discrepancy  could  be  explained  in  one  of  two  ways: 

(a)  the  amplitudes  of  the  two  fields  are  similar,  but  they  are  out  of  phase 
with  each  other,  or  (b)  the  corrected  fields  are  more  like  the  Flattery 
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analysis  fields  than  they  are  like  the  observations,  so  that  comparing  the 
Flattery  analyses  with  the  corrected  fields  is  not  like  comparing  them  with 
the  observations.  Looking  back  at  Table  1,  there  is  certainly  some  evidence 
that  most  of  the  corrected  fields  have  only  been  brought  less  than  half  way 
towards  the  observations  from  the  Flattery  analyses.  Thus,  they  are  for  the 
most  part  more  a  product  of  the  Flattery  analysis  than  they  are  of  the  obser¬ 
vations.  However,  in  order  to  examine  the  possibility  of  a  phase  error,  the 
absolute  value  of  the  error  between  the  corrected  fields  and  the  Flattery 
analyses  must  be  computed.  Values  for  the  globally  averaged  absolute  value 
of  the  difference  between  the  two  fields  are  given  in  Table  8,  along  with 
the  differences  between  the  global  averages  of  corrected  fields  and  Flattery 
fields.  The  units  for  the  quantities  in  the  first  two  columns  are  given, 
and  the  values  in  parentheses  are  the  ratio  of  the  magnitude  of  the  values 
in  the  second  column  to  those  in  the  first.  Note  that  the  absolute  error 
between  the  two  fields  is  small  compared  to  the  average  magnitude  of  each 
quantity  (column  4)  especially  for  T,  Z,  and  q,  so  that  a  relatively  small 
phase  shift  occurs  in  the  corrected  field  with  respect  to  the  Flattery  field. 
However,  notice  that  the  bias  ([c]  -  CF ] )  values,  though  they  are  small  in 
magnitude,  make  up  a  significant  fraction  of  the  total  magnitude  difference 
(column  1)  between  the  fields  for  those  quantities  (T,  Z,  and  q)  which 
showed  the  most  significant  amount  of  correction  as  seen  in  Table  1.  It  is 
concluded  from  these  results  that  changes  in  the  Flattery  fields  due  to  the 
successive  corrections  were  relatively  small,  but  were  appreciably  effective 
in  biasing  the  fields  towards  the  observations,  especially  for  T,  Z,  and  q. 

TABLE  8.  GLOBAL  AVERAGE  OF  DIFFERENCES  BETWEEN  CORRECTED  AND  FLATTERY  FIELDS 
FOR  JANUARY  WEEK 


850  u 

v 
T 
Z 

q 

500  u 

v 

T 

Z 

q 


[lci  - 

FJ1 

Cc]  -  [F] 

[C]  -  [F] 

[ICf  -  F 

[|C.  -  F.|] 

[|F|] 

.909 

-1 

ms_i 

-.165 

(.18) 

.141 

.855 

ms 

-.036 

(.04) 

.242 

.537 

°K 

-.247 

(.46) 

.002 

4.192 

9P«> 

.707 

(.17) 

.003 

.387 

g  kg 

-.202 

(.52) 

.056 

.924 

-.006 

(.01) 

.084 

.896 

.040 

(.05) 

.104 

.315 

-.166 

(.53) 

.001 

6.223 

2.532 

(.41) 

.001 

.118 

-.044 

(.38) 

.065 

refgga?*  "zw.ttr.#nz mtm 
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In  summarizing  the  comparison  between  the  Flattery  fields  and  the  cor¬ 
rected  fields,  it  is  evident  that  the  corrected  fields  still  have  many  of 
the  physical  and  spectral  characteristics  of  the  Flattery  fields.  Both  the 
absolute  difference  averaged  over  grid  points  and  the  difference  in  appor¬ 
tionment  of  power  were  small  between  the  two  fields  for  all  of  the  quanti¬ 
ties.  That  the  phase  error  is  small  is  a  tribute  to  the  Flattery  scheme  in 
that  highs  and  lows  and  minima  and  maxima  had  the  same  locations  as  indi¬ 
cated  by  the  observations.  Furthermore,  the  relative  amplitudes  of  the 
waves  were  generally  preserved  by  the  analysis.  The  corrections  acted  in 
such  a  way  as  to  draw  the  average  value  of  the  field  more  toward  that  of  the 
observations  (especially  for  T,  Z,  and  q)  rather  than  change  the  fields  at 
particular  locations. 

V.  Summary  and  Conclusions 

The  following  is  a  summary  of  the  similarities  and  differences  noted 
when  comparing  the  spectral  characteristics  of  global  specific  humidity  and 
layer  precipitable  water  fields  with  those  of  zonal  and  meridional  wind  com¬ 
ponents,  temperature,  and  geopotential  height: 

a.  The  variation  of  the  zonal  mean  power  of  T,  Z,  q,  W  with  latitude  is 
similar  for  all  four  quantities,  with  a  maximum  at  the  equator  and  a  mini¬ 
mum  at  the  poles.  The  zonal  mean  for  q  and  W  is  equal  to  the  global  mean 

at  latitudes  lower  than  for  T  and  Z,  indicating  a  larger  meridional  gradient 
of  the  zonal  mean  for  the  moisture  quantities.  The  latitudinal  variation  of 
the  zonal  mean  of  u  reflects  the  distribution  of  easterlies  and  westerlies 
and  thus  oscillates  several  times  between  the  poles. 

b.  The  latitudinal  distribution  of  total  zonal  eddy  power  reveals  a  bimodal 
character  for  T  and  Z,  with  maxima  at  mid- latitudes  and  minima  at  the  equa¬ 
tor  and  both  poles.  This  is  also  true  for  u  and  v,  although  to  a  somewhat 
less  pronounced  extent.  The  moisture  variables  have  a  zonal  eddy  power 
which  is  a  broad  maximum  at  equatorial  latitudes  and  a  minimum  at  both  poles. 
It  appears  that  this  latitudinal  structure  is  primarily  due  to  the  satura¬ 
tion  specific  humidity  acting  as  a  more  severe  limit  to  the  variation  about 
the  zonal  mean  in  colder  regions  than  it  is  in  warmer  regions  because  of  its 
lower  value  at  lower  temperatures . 

c.  The  zonal  harmonic  analysis  results  showed  that  the  power  was  more  re¬ 
stricted  to  lower  wave  numbers  for  u,  T,  and  Z,  and  more  spectral  spreading 
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was  evident  in  v,  q,  and  W.  The  T  and  Z  spectral  distribution  varied  more 
significantly  with  latitude,  with  a  greater  restriction  of  power  to  lower 
wave  numbers  at  mid- latitudes  and  more  spectral  spreading  in  the  tropics. 
Power  shifts  to  lower  wave  numbers  with  height  were  evident  in  the  January 
case  for  q  and  W  and  in  the  July  case  for  T  and  Z. 

d.  In  the  spherical  harmonic  results,  the  wind  components  show  the  great¬ 
est  spectral  spreading,  with  most  of  their  power  in  the  composite  wave  num¬ 
ber  3-10  band.  They  also  have  a  significant  amount  of  their  power  in  the 
11-20  wave  number  band.  This  is  unlike  the  power  distribution  for  the  other 
four  variables,  which  have  a  much  smaller  amount  of  power  (less  than  10  per¬ 
cent)  contributed  by  wave  numbers  greater  than  10.  There  is  very  little 
variation  in  power  distribution  with  height  in  u  and  v.  The  q  and  W  spec¬ 
tra,  on  the  other  hand,  show  power  spreading  to  higher  wave  numbers  with 
height,  while  power  is  more  concentrated  at  lower  wave  numbers  with  height 
for  Z  and  to  a  lesser  extent  for  T.  The  power  for  T  results  primarily  from 
the  number  2  meridional  wave  of  the  zonal  mean,  and  there  is  very  little 
change  in  this  with  season  or  height.  In  general,  T  and  Z  have  most  of  the 
power  in  the  lowest  wave  numbers  (0-2) ,  q  and  W  have  roughly  an  equal  amount 
of  power  contributed  from  low  and  medium  (3-10)  wave  numbers,  and  u  and  v 
have  most  of  their  power  from  the  medium  wave  numbers,  with  a  significant 
contribution  from  the  high  (11-20)  wave  numbers. 

e 

e.  while  the  meridional  wave  number  2  component  of  the  zonal  mean  is  the 
dominant  spherical  component  for  T,  Z,  q,  and  W,  the  results  showed  that  the 
meridional  wave  number  1  component  of  the  zonal  mean  grew  at  the  expense  of 
wave  number  2  in  the  July  case  for  T  and  z.  Thus,  the  zonal  mean  values  in 
the  Northern  Hemisphere  in  July  were  higher  than  their  counterparts  in  the 
Southern  Hemisphere  in  January,  making  the  July  zonal  mean  temperature  plot 
much  more  asymmetric  about  the  equator  than  the  January  plot.  The  higher 
Northern  Hemisphere  summer  temperatures  were  attributed  to  the  larger  land 
to  ocean  ratio  in  that  hemisphere.  This  asymmetry  lowered  the  meridional 
wave  number  2  contribution  from  q  and  W  as  well. 

The  results  presented  in  this  Btudy,  to  the  extent  that  they  are  con¬ 
sidered  an  accurate  portrayal  of  the  spectral  characteristics  of  the  global 
moisture  distribution,  are  intended  to  be  used  as  a  baseline  for  testing  the 
reliability  of  an  objective  analysis  scheme  for  moisture.  Clearly,  the 
unique  characteristics  of  moisture  must  be  faithfully  reproduced  in  the 
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gridded  analysis  £ield,  and  not  be  distorted  by  the  analysis  method  itself. 
This  implies  the  need  for  modification  of  a  mass-motion  objective  analysis 
scheme  to  make  it  useful  for  moisture  analysis,  or  development  of  an  entirely 
separate  scheme  for  moisture  analysis. 
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